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We present detailed measurements of the longitudinal resistivity pxx{T, H), and the upper critical 
field Hc2 of Nd(Oo.7Fo.3)FeAs single crystals in strong DC and pulsed magnetic fields up to 45 T 
and 60 T, respectively We found that the field scale of Hc2 is comparable to Hc2 100 T of 
high Tc cuprates. Hc2{T) parallel to the c-axis exhibits a pronounced upward curvature similar to 
what was extracted from earlier measurements on poly crystalline La(0,F)FeAs, Nd(0,F)FeAs and 
Sm(0,F)FeAs samples. Thus this behavior of H^{T) is indeed an intrinsic feature of oxypnictides, 
rather than manifestation of vortex lattice melting or granularity. The orientational dependence 
of Hc2{0) as a function of the angle between H and the c-axis shows deviations from the one- 
band Ginzburg-Landau scaling. The mass anisotropy parameter 7(T) = {mc/rriahY^'^ — HI2/HC2 
obtained from these measurements decreases as temperature decreases, from 7 ~ 9.2 at 44 K to 7 ~ 5 
at 34 K, where || and ± correspond to H parallel and perpendicular to the ab planes, respectively. 
Spin dependent magnetoresistance and nonlinearities in the Hall coefficient suggest contribution 
to the conductivity from electron-electron interactions modified by disorder reminiscent that of 
diluted magnetic semiconductors. The Ohmic resistivity pxx{T, H) measured below Tc but above the 
irreversibility field exhibits a clear Arrhenius thermally activated behavior p = po exp[— ^a(T, H)/T] 
over 4 — 5 decades of pxx- The activation energy H/a 

(T, H) has very different field dependencies for 
H\\ab and H ± ab, varying from 4x 10^ K at if = 0.2 T to - 200 K at if = 35 T. We discuss to what 
extent different pairing scenarios suggested in the literature can manifest themselves in the observed 
behavior of iic2, using the two-band model of superconductivity in oxypnictides. The results indicate 
the importance of paramagnetic effects on Hc2{T) in oxypnictides, which may significantly reduce 
iic2(0) as compared to iic2(0) ^ 200 — 300 T based on extrapolations of Hc2{T) near Tc down to 
low temperatures. 

PACS numbers: 74.70.-b,74.25.Fy,74.72.-h,74.81.Bd 



The recently discovered layered superconducting oxyp- 
nictides with high transition temperatures T^^ are based 
on alternating structures of FeAs layers and REG lay- 
ers of rare earth (RE) oxides. Similar to the high- 
temperature superconducting cuprates, superconductiv- 
ity in oxypnictides emerges upon doping of a parent an- 
tiferromagnetic state. As the REG planes are doped, the 
ionically bonded REG donates an electron to the cova- 
lently bonded FeAs plane^, suppressing the global anti- 
ferromagnetism and resulting in superconductivity. 

Previous low ^eld^d,^^ and high field transport 
studies^'^ on different polycrystalline oxypnictides have 
shown, that their magneto-transport behavior is inter- 
mediate between conventional low- Tc superconductors 
and the high- temperature cuprates. However, while 
La(G,F)FeAs behaves as an intermediate-Tc supercon- 
ductor similar to MgB2 in which thermal fluctuations of 
vortices do not signiflcantly affect the H-T diagram to 
the extent that they do in the layered cuprates, the higher 
Tc oxypnictides like SmFeAsGo.ss and Nd(G,F)FeAs ex- 
hibit larger mass anisotropics, enhanced thermal fluctu- 
ations, and a Ginzburg parameter comparable to that 



of YBGG. Thus, the question arises, whether this dis- 
tinct behavior stems from different critical tempera- 
tures, the presence of magnetic ions^ii^iii, or is due to 
the electromagnetic granularity of these polycrystalline 
samples^^^^^^^^. For example, is the pronounced upward 
curvature of Hc2{T) along the c-axis extracted from mea- 
surements of the resistivity onsets of polycrystals^ in- 
deed an intrinsic feature of H^{T) in oxypnicties, or 
does it reflect an extrinsic effect of magnetic granular- 
ity, the onset of the irreversible critical state behavior 
and/or melting of the vortex lattice? Moreover, given 
that the nature of superconductivity in the oxypnic- 
tides is still poorly understood, several different sce- 
narios have been proposed. In particular, many the- 
oretical models are based on the multiband electronic 
structure of the oxypnictides and the fact that the su- 
perconducting state emerges upon doping of the par- 
ent antiferromagnetic semimetaU^ii^iiiii^ii^. The multi- 
band effects would manifest themselves in the distinct 
temperature and orientational dependencies of Hc2{T), 
as has been well documented in detailed studies of 
MgBs^-^^^^^^^^^^^^^. To address these 
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issues, we made detailed measurements of the magneto- 
transport properties of a Nd(Oo.7Fo.3)FeAs single crystal 
in very high magnetic fields. 



I. SAMPLES 

The crystal was made by the flux method using NaCl 
as the flux. First the starting Nd (purity 99.95 %) and 
As (purity 99.99 %) were mixed in 1:1 ratio, ground and 
pressed into a pellet shape. Then it was sealed in an evac- 
uated quartz tube and reacted at 800 °C for 10 hours. 
The resultant pellet was ground together with a stoi- 
chiometric powder mixture of NdFa (purity 99.95 %), 
Fe203 (purity 99.9 %) and Fe (purity 99.9%) resulting in 
Nd(Oo.7Fo.3)FeAs compound. Again it was pressed into 
a pellet and put together with NaCl powder of mass ratio 
(NaCl : Nd(0,F)FeAs = 10:1) and sealed in an evacuated 
quartz tube and reacted at about 1050 °C for 10 days. 
Then it was cooled down at a rate of 3 °C/hour to 850 
°C, followed by a quick cooling to room temperature by 
shutting off the power of the furnace. The resulting prod- 
uct contains mainly plate-like small crystals with lateral 
dimensions of 5-50 /im and ~ 1 /am thickness. 

Previous extensive transport^^^^^'^^ and heat 
capacity^ studies on Nd(Oo.82Fo.i8)FeAs single 
crystal revealed the onset of superconductivity at 
Tc ^ 47 K, the upper critical fields H^^iO) = 304 T and 
H^{0) = 62 — 70 T extrapolated from the Werthamer- 
Helfand-Hohenberg (WHH) formula^^. These samples 
also have a relatively small mass anisotropy parameter 
7 = {rUc/maY^'^ ^ 5 where rua and rric are the effective 
masses along the ab plane and the c-axis, respectively. 
However, the above studies were performed in magnetic 
fields not exceeding 9 T, where most of the multi-band 
effects would not manifest themselves in the Hc2{T) 
curves. 



II. RESISTIVE TRANSITIONS AND UPPER 
CRITICAL FIELDS 

For our experiment we used three different high-field 
magnets at the National High Magnetic Field Labora- 
tory (NHMFL): 35 T DC resistive and 45 T hybrid mag- 
nets at Florida State University and 65 T pulsed field 
magnet at Los Alamos National Laboratory. The inset 
to Fig. [T](a) shows the scanning electron micrograph of 
the Nd(0,F)FeAs sample with six Pt contacts fabricated 
by the focused ion beam (FIB) technique. The distance 
between voltage contacts is ~ 13 /im, and the sample 
thickness is ~ 1.4 jam. The sample was glued to a conduc- 
tive silicon substrate to ensure good thermalization dur- 
ing measurement, and to prevent electric charge build-up 
during the FIB procedure. The sample was placed on a 
rotating platform, which allows us to change its orienta- 
tion with respect to the magnetic field in situ. The longi- 
tudinal resistivity pxx and the Hall coefficient Rh in high 
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FIG. 1: Longitudinal resistivity pxxiT) at if = and H — 
35 T for H\\ah and H ± ab. At B = superconductivity 
onset is observed at Tc :^ 46 K. Inset shows scanning electron 
micrograph of the sample. 



magnetic fields were measured using a lock-in technique 
with AC excitation current below 0.3 mA at frequency 22 
Hz. Pulsed field measurements were conducted using in- 
house synchronous digital lock-in technique at frequency 
of 16 kHz and similar levels of current. 

Figure [H shows pxxiT) as a function of temperature 
T at if = and 35 T for different sample orientations. 
The magnetic field perpendicular to the ab plane shifts 
the superconducting transitions to lower T, reducing Tc 
(from 46 K to 21 K at if = 35 T) without significant 
change of the shape of pxx(T). This behavior is reminis- 
cent of magnetotransport in low-T superconductors. In 
contrast, fields applied parallel to the ab plane change 
the shape of the resistive transitions more than they do 
for H\\c^ but less than what was observed on polycrys- 
talline Nd(0,F)FeAs^. Figure [1] shows that the pxx{T) 
dependence is sublinear at higher T but superlinear below 
T < 100 K. Over the same temperature range, the Hall 
resistivity pxyiT) shown in Fig. [2] has a relatively strong 
temperature dependence for such an optimally doped ma- 
terial. This dependence is much stronger than was ob- 
served on a polycrystalline sample with much lower dop- 
ing, X = 6 with pxy{T) exhibiting a kink around 
T ^ 100 K. Such temperature dependencies of the trans- 
port coefficients near ~ 100 K may result from a multi- 
band conductivity or a structural transition. The inset 
to Fig. [2] shows pxy{H) at selected temperatures. On this 
field scale, pxy{H) does not exhibit field nonlinear it ies, 
at 40 < T < 300K, similarly to polycrystalline samples^. 
However, a strong dependence of Rh{H) is observed in 
this single crystal at lower temperatures as will be dis- 
cussed later. 

Figures [3]-[6] show the detailed field and temperature 
dependencies of pxx{T^H)^ from which the upper crit- 
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FIG. 2: Hall coefficient Rh as a function of T measured in 
superconducting magnet at ±9 T (solid line), and determined 
from linear fits to pxy(H) for —35 < if < 35 T or —9 < 
H < 9 T (solid squares), as shown in the inset. Rh for 
poly crystalline Nd(Oo.94Fo.o6)FeA sample is also shown. 




FIG. 3: pxx as a function of temperature at various magnetic 
fields for H\\ab. 
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FIG. 4: pxx as a function of temperature at various magnetic 
fields for H A_ ab. Dotted lines show how the onset of the 
transition was determined for H — 20 T. Excitation current 
/ < 0.3 mA, at frequency / = 22 Hz was applied in the ab 
plane. 
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FIG. 5: The longitudinal resistivity pxx{H) as a function of 
magnetic field in for H\\ab measured in pulsed magnet up to 
60 T. 



ical fields parallel and perpendicular to the ab planes 
were extracted. Here Hc2 for both field orientations was 
measured at the onset of the superconductivity transi- 
tion defined by the intersection point from linear ex- 
trapolations of pxxiT^H) T < Tc and pn{T^H) at 
T > Tc, as shown in Fig. H Here Pn{T,H) is 
the resistivity in the normal state. We also extracted 
the mid transition field H^o{T) and the resistivity on- 
set fields Hio{T) and Ho.5{T) defined by p^^{T,H) = 
0.bpn{T,H), p,,{T,H) = 0.1pn{T,H), and p,,{T,H) = 
0.05pn{T^ H)^ respectively. The fields Ho_^{T) are close 



to the irreversibility field if*, which quantifies the onset 
of vortex critical state. The resulting temperature de- 
pendencies of all these fields are summarized in Figs. [8] 
and[9l Notice that the R{T) curves shown in Figs. [3]- 
[6] exhibit a small peak just before the transition which 
disappears at > 4T for H||c and H > 6T for H||a6. 
Such resistance peak anomaly has been observed in other 
superconductors and is usually ascribed to the effect of 
paramagnetic ions, sample inhomogeneities, contact ar- 
rangements, etc.^-— i^i^. Since the mechanism of this 




FIG. 6: The longitudinal resistivity pxx{H) as a function of 
magnetic field in for H J_ ab measured in DC resistive magnet 
up to 35 T at various temperatures. 
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FIG. 7: Magnetoresistance measured at T = 50 K for two 
different sample configurations. Squares denote data from 
temperature sweeps at fixed magnetic field, while solid lines 
were measured at fixed T. (b) The same data plotted versus 
H^. Dotted lines show linear fits taken for if ^ 15 T. (c) 
Data plotted after subtraction of the contribution, (d) 
Hall coefficient Rh{H) = pxy{H)/H versus magnetic field 
measured at T = 65 K . 
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FIG. 8: Temperature dependence of the upper critical field 
i7j,2(T) along the ab plane obtained from dc (filed symbols) 
and pulsed (open symbols) measurements. The data ex- 
tracted from the results shown in Fig. O and Fig. O show 
the temperatures at which the resistance reaches 0.5 %, 10 %, 
50 % of the normal state resistance, as extrapolated linearly 
from the pn{T,H) temperature dependence above Tc{H). 
The onsets of superconducting transition were determined as 
shown in Fig. |4] and correspond well to the 90 % transitions. 
The dotted lines are guides to the eye. 
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anomaly is still not completely understood, we did not 
take it into account when extracting Hc2 from the data 
shown in Figs. [3]-[6l 

Figure [TJa) shows the magnetoresistance (MR) at 
T = 50 K, just above the onset of superconductivity. 
The MR is positive for i^||c and negative for for i^||ab. 
The difference could be attributed to orbital origin of 
this MR. As shown in Fig. [T^b), the magnetoresistance, 



FIG. 9: Temperature dependence of the upper critical fields 
H^(T) along the c-axis obtained from the dc measurements. 
The data extracted from the results shown in Fig. |4l and 
Fig. [51 show the temperatures at which the resistance reaches 
0.5 %, 10 %, 50 % of the normal state resistance, as extrap- 
olated linearly from the pn{T,H) temperature dependence 
above Tc{H). The onsets of superconducting transition were 
determined as shown in Fig.|4]and correspond well to the 90 % 
transitions. The dashed lines are guides to the eye. 
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{pxx{H) — pxx{^)) / Pxx{^) oc H'^ for both field orienta- 
tions and i7 ^ 15 T exhibits a quadratic field depen- 
dence. However, if this quadratic contributions are re- 
moved, a weak remanent positive MR is observed at low 
fields ^ 15 T, as shown in Fig. [Tl^c). This remanent 
MR is identical for both field orientations, which in turn 
suggests a spin mechanism behind this positive MR com- 
ponent. Interestingly, such a spin dependent positive MR 
is ubiquitous in diluted magnetic semiconductors (DMS) 
in the paramagnetic phase^^^^^'^^. It originates from the 
giant spin-splitting of the electron states, which con- 
siderably affects quantum corrections to the conductivity 
brought about by the effect of disorder modified electron- 
electron interactions^^. Since in DMSs is proportional 
to the magnetization M of the localized spins, this pos- 
itive MR scales with B and T like the Brillouin func- 
tion. Importantly, this positive MR is absent in a fer- 
romagnetic phase, when localized moments are already 
aligned at = 0. This mechanism can also produce 
MR in Nd(0,F)FeAs, provided that at least some of lo- 
calized magnetic moments ji (most likely of Nd^+ ions 
with ji = 3.6/iB) are in a paramagnetic phase. At the 
same time, as shown in Fig. [Tl^d), the Hall coefficient 
Rh strongly depends on H at low field. Such a depen- 
dence could also result from quantum corrections to pxy 
from electron-electron interactions, but maybe a result of 
multiband conductivity, as well. 

Based on the results shown in Figs. [8] and [9l the fol- 
lowing points can be made: 

1. The 100-200 T field scale of Hc2 in 
Nd(Oo.7Fo.3)FeAs is comparable to that of high-Tc 
cuprates. 

2. H^{T) perpendicular to the ab plane exhibits a 
pronounced upward curvature similar to what was first 
reported for polycrystalline La(0,F)FeAsi. This is very 
different from the one-band WHH behavior and appears 
to be an intrinsic property of oxypnictides rather than 
manifestations of the vortex melting or the onset of pin- 
ning at = H%T). 

3. The slopes H'JI = \dH^}^/dT\ and H'^^ = \dH^^/dT\ 
for both field orientations increase significantly just a few 
Kelvins below Tc. Thus, low- field < 9 T measure- 
ments may underestimate the actual values of — 
^o/^^il^Tc and H% = ^o/^^U^coTc at T ^ T„ from 
which the coherence lengths ^^(T) = (1 - T/Tc)-^^^^aO 
and £,c{T) = (1 — T/Tc)~^/^^cO in the ab plane and along 
the c-axis can be extracted. As it is evident from Figs. [8] 
and[9l the high-temperature slopes ^ci attain 
different values above :^ 10 T, so the use of high mag- 
netic fields becomes decisive for revealing the true in- 
trinsic superconducting parameters of Nd(Oo.7i^o.3)FeAs 
single crystals. This may also explain why low-field 
measurement^i2ii2^ show lower values of dHl2/dT and 
dH^2/dT. 



III. THERMALLY ACTIVATED RESISTIVITY 

Shown in Figs [10] and [11] are the temperature de- 
pendencies of the ohmic resistivity pxx{T^H) measured 
above the irreversibility field, i^* < H < Hc2 for both 
field orientations. These data can be described well by 
the Arrhenius dependence 

Inp,, = lnpo-^a(T,i/)/T, (1) 

characteristic of thermally-activated flux flow (TAFF) of 
vortices, which has been studied extensively on high-Tc 
cuprates. Here Ea{T^B) = —d\npxx/d{l/T) is the acti- 
vation energy determined by hopping of vortex bundles 
in the pinning potential. Several features of the data 
shown in Figs. [TO] and [TT] should be pointed out. 

1. The Arrhenius dependence holds over 4-5 decades 
of pxx^ down to the limit of sensitivity of our transport 
measurements. This indicates that thermally-activated 
vortex hopping in Nd(Oo.7Fo.3)FeAs is similar to that of 
high-Tc cuprates and is much more pronounced than in 
lower- Tc pnictides like La(Oi_ajFa^)FeAs. The good lin- 
ear behavior evident from Figs [10] and [11] also indicates 
that the low-temperature dependence of Ea (T) is approx- 
imately linear, EaiT) = Em{^ — T/Tm)^ where Tm is a 
parameter with the dimensionality of temperature. In 
this case the term T/Tjn can be eliminated by re-defining 
Po po^'^'p{Em/Tm) in Eq. ([T|) without changing any 
observed characteristics. 

2. The higher temperature log pxx{T) data level off 
at a field independent value corresponding to the normal 
state resistivity pn- 

3. The \ogpxx{T) lines for different H\\c extrapolate 
to the same temperature, ^ Tc, indicating that T^ = Tc, 
Po ^ Pn exp(£^^Tc) and Em is practically independent of 
H. However, for the parallel field orientation, this is not 
the case because of the field dependence of Em, as shown 
below. 

From the low-temperature slope of logp^^^^ the 
field dependencies of the activation energy Ea = 
—d\npxx/d{l/T) can be obtained. The results pre- 
sented in Fig. [12] for both field orientations show that 
Ea (3 - 4) X 10^ K at low fields. The acti- 
vation energy Ea — 2 x 10^ K for a polycryst aline 
Nd(0,F)FeAs measured previously^ is higher than for the 
single crystals studied in this work (similar thermally- 
activated resistivity was recently observed on polycrys- 
talline Nd(Oo.82Fo.i8)FeAs«^). The field dependencies of 
Ea{B) for i^||c and H ±c are also markedly different. 
Ea{H) for i^||c exhibits a rather conventional field de- 
pendence characteristic of thermally- activated transport: 
a constant value at low fields < 3 T where single- 
vortex pinning dominates, followed by a power law de- 
crease Em oc H~^-^ characteristic of collective creep at 
higher fields, H > 3 T— . By contrast, the activation en- 
ergy for H\\ab exhibits a rather weak power law decrease 
Em cx H~^'^^ in the entire field interval where pxx{T, H) 
was measured. 
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FIG. 10: Arrhenius plot for pxxiT) at different magnetic fields 
perpendicular to the ab plane. 
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FIG. 12: Field dependencies of the activation energy Ea{H) 
for H J_ ab and H\\ab. The solid squares show the data for 
Nd(0,F)FeAs poly crystal measured earlier. 
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FIG. 11: Arrhenius plot for pxx{T) at different magnetic fields 
parallel to the ab plane. 
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FIG. 13: The resistivities pxx{T,H) measured at 44 K for 
different angles between H and the ab plane. 



The field and temperature dependencies of Ea{T^H) 
can be described by the following scaling function 



IV. ANGULAR DEPENDENCIES 



Ea 



E^(l-T/T,Y 
[l^H/Ho{T)f 



H 



(2) 



where Hq quantifies the field above which vortex dy- 
namics is determined by collective effects of vortex 
interaction^^. From the data presented above, we ob- 



tain Hq 3 a 
for mab. 



1, and f3 ^ 1.1 for H||c and (3 ^ 0.17 



The TAFF resistivity pxx{T, H^O) was measured at 
different angles between H and the ab plane, as il- 
lustrated by a representative example shown in Fig. [131 
The angular dependencies of Hc2{0) extracted from the 
PxxiT, H,e) at 90 %, 50 %, and 10 % of Pn{T,H) are 
shown in Fig. [TH 

It is clearly seen that the anisotropy as 7 = h12/H^ 
decreases when calculated at lower resistive transitions. 
At the same time i^^2^''(^) shown in Fig. [TH notably de- 
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FIG. 14: Angular dependence of Hc2{0) defined at pxx — 
0.9 0.5 O.lpn at 44 K. The solid lines show the fits with Eq. 



partures from the standard Ginzburg-Landau formula: 



c2 



Vsin^ 



7 ^ cos^ I 



(3) 



where 7 = {mc/rriaY^'^ is defined by the ratio of effec- 
tive masses along the ab plane and the c-axis. The fit 
to Eq. ([3|) yields 7 ~ 8.9 in this case. Data taken at 
50 % and 10 % departure from Eq. ([3]) even stronger, 
which may result from increasing contribution of TAFF 
transport at lower resistivities. 

The temperature dependencies of 7(T) obtained from 
the Hc2 data shown in Figs. [8] and [9] for iYio, i^so and Hc2 
are shown in Fig. 15. Contrary to the standard one-band 
behavior, 7 turns out to be temperature dependent, the 
values of 710 (T), 750 (^) and 7(T) being rather different, 
especially just below Tc. While the parameter of Hc2 
anisotropy, 7(T) = H^^2/ ^^2 decreases from ^ ^ dX Tc to 

5 at 34 K, the field Hiq becomes more anisotropic upon 
decreasing T. Below T ^i:^ 38 K, however, these values of 
7(T) converge. 



V. Hc2 IN TWO-BAND SUPERCONDUCTORS 

To understand how different pairing scenarios could 
manifest themselves in the observed temperature depen- 
dencies of Hc2{T), we use the two-band BCS theory in 
which Tc for interband scattering is given by^i^ 



Tco l^exp[-(A+ - sXo)/2w], 



(4) 



where A± = An ± A22 , ^ = An A22 - A12 A21 , 5 = sign(i(;) , 
and Ao = (A?_ + 4Ai2A2i)^^^, A12 and A22 are pairing 
constants in bands 1 and 2, and A12 and A21 quantify 
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FIG. 15: Temperature dependencies of 710, 750 and 790. 
The open squares show the data from Ref. 35 measured on 
Nd(Oo.82Fo.i8)FeAs single crystal with higher Tc — 51.5 K. 
These data are plotted at reduced temperatures. 



interband coupling. For s-wave pairing, intraband impu- 
rity scattering does not affect Tco, while weak interband 
scattering reduces Tc according to ^^i^^i^^ 



Tr = T 



cO 



8A0 



[A07+ + -^-7- - 2A21712 - 2A12721] (5) 



where 7± = 712 ± 721, and 712 and 721 are scatter- 
ing rates between bands 1 and 2, A/'i7i2 = N2J21 and 
^1^12 = ^^2^21, where Ni and A^2 are partial densities 
of states in bands 1 and 2, respectively. In the following 
we discuss two scenarios: 1. Conventional s-wave pair- 
ing dominated by strong intraband coupling, w > 0. 2. 
Strong interband pairing w < which can result in the 
TT shift between the order parameters on two bands for 
interband repulsion A12 < (s^ pairing^^). 

As follows from Eq. ([5|), interband repulsion with 
A12 < and A21 < enhances Tc depression due to inter- 
band impurity scattering between different pieces of the 
Fermi surface with opposite signs of the order parame- 
ter. This suppression may be reduced by strong coupling 
effects^^. By contrast, for strong interband attraction, 
4A12A21 > A^, Eq. ([5D yields 



7r7i2 




(6) 



Here Tc suppression is much 
pairing because the factor (1 



weaker than for the 
^JNl/N2f reduces the 
effect of interband scattering, making Tc independent 
of 712 in the special case of A^i = A^2- For instance, 
for Ni = Eq. ([6|) gives a rather weak suppres- 

sion of Tc — Tco(l — 0.02^) even in the dirty limit, 
g = 7i2/2^T,o > 1. 
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In the dirty limit, the equation for Hc2^ which takes 
into account both orbital and Zeeman pairbreaking for 
negligible interband scattering, can be written in the fol- 
lowing parametric form^: 

\iit = -[Ui{h) + U2{h) + Xo/w]/2 + (7) 
s[{Uiih) - U2{h) - X-fwf/A + Ai2A2i/^Y/', 

Hc2 = 2(l)okBT,th/hDo, (8) 

where t = T/Tco, and the functions Ui^2 are defined by 

Ui,2{h) = Re^[l/2 + (i + D,^2/Do)h] - V^(l/2), (9) 

is the di-gamma function, Di and D2 are diffusiv- 
ities in band 1 and 2, I^o = ^/2m, and the parameter 
h runs from to 00 as T varies from Tc to 0. For equal 
diffusivities, rj = D2/ Di = 1^ and negligible Zeeman pair- 
breaking, 1^0 <^ ^1,2, Eq. ([8]) simplifies to the one-band 
de-Gennes-Maki equation, Int -\- U{h) = 0. If the mag- 
netic field H is inclined by the angle with respect to the 
ab planes, Eqs. ([8])- ([9]) also describe the angular depen- 
dence of Hc2{0) if the intraband diffusivities are replaced 
by their angular-dependent values: 

D,r^{0) = [D^^^^ sin^ + D^^^Di^^ cos^ 0]'^^ (10) 

where m = 1,2, and Din' and are the in-plane and 
the c-axis principal values of . 

Hc2{T) curves shown in Fig. [8] and [9] exhibit pro- 
nounced upward curvature for the field along the c-axis, 
and the downward curvature for the field along the ab 
plane, similar to the behavior of Hc2{T) first observed 
on polycrystalline LaOo.sgFo.iiFeAs^. Such temperature 
dependence of Hc2{T) has been often observed on MgB2, 
suggesting that the two-band theory in which two bands 
have different diffusivities can also be used to describe 
Hc2{T) in oxypnictides (the superconducting gaps on the 
disconnected pieces of the Fermi surface can be either 
different or the same). If interband scattering is negli- 
gible, neither Tc nor Hc2 depend on the sign of A12, so 
we consider two different cases: 1. dominant intraband 
coupling, w = A11A22 — A12A21 > and 2. dominant 
interband coupling, w = A11A22 — A12A21 < 0. 

Shown in Fig. [T6lf a) is an example of the fit of Eq. ([8])- 
(|9|) to the data for case 1 with An = A22 = 0.5, and A12 = 
A21 = 0.25, in which case Tc = 1^ exp[— l/(Aii + A12)] = 
50K corresponds to 1] 190K. Fig. [THlf a) shows the 
fit without taking paramagnetic effects into account for 
D2 = 0.081^1, which suggests that scattering in one band 
is significantly stronger than in the other. We do not dis- 
cuss here microscopic mechanisms behind such band dis- 
parity; instead we focus on the set of parameters which 
fit the data and the extent to which these parameters are 
different for the case of w > and w < 0. Extrapola- 
tions of i^j2(^) to T = based on the observed slope at 
Tc suggests H^c2i^) ^ 220T, well above the weak coupling 
BCS paramagnetic limit Hp[T] = 1.84Tc[K] 90T. 

Fig. [T6lfa) shows the fit in the case of strong interband 
pairing A12A21 > A11A21 for An = A22 = 0.49, A12 = 



A21 = 0.5 and D2 = O.OOTD^. Reproducing the observed 
upward curvature of (T) for this case requires a much 
greater disparity of the band diffusivities than for the 
previous case. Reducing An and A22 shifts the region 
of upward curvature to lower temperatures, which makes 
fitting the data even more difficult. These features may 
put serious constraints on the pairing scenario. 

Given the very high extrapolated values of i^j2(0) 
in Fig. [T6lf a). the role of paramagnetic effects becomes 
very important. Shown in Fig. [T6lfc) is the fit for 
the same parameters as in Fig. flUT a) but with the ac- 
count of the Zeeman pairbreaking for d[^^ = 10 Dq and 

[L)^''^L>^'V^i''^^i'^^]^^^ = 0-7^0. In this case the values 
of Hc2{T) at low temperatures are significantly reduced 
as compared to what may be expected from the fit in 
Fig. [T6lf a). which ignores the paramagnetic effects. 

Several remarks should be made regarding the fits in 
Fig. [T6lf a-c). First, applying weak coupling Eqs. ([8|)-([9|) 
for a superconductors with Tc = 50 K obviously ne- 
glects the strong coupling renormalization effects of the 
Eliashberg theory. However given the lack of microscopic 
theory of superconductivity in oxypnictides, the use of 
the BCS-type approach may have its own merits given 
that the shape of the Hc2{T) curve described by Eqs. 
(|8|)-(|9]) depends only on the ratio of the constants Xmn 
but it is basically insensitive to the particular s-wave 
pairing mechanisms. The usual way of treating strong 
coupling effects in Eqs. ([8])- ([9]) is to express them in 
terms of experimentally observed parameters^ such as 
Tc and Di and D2 extracted from normal state trans- 
port measurements and the observed slopes dHc2/dT. 
The second remark is that the actual paramagnetic limit 
Hp (1 X)Hp^^ is significantly enhanced by strong 
coupling effects as compared to the BCS value^^. Here 
A ~ 1 is a characteristic bare coupling constant for the 
relevant exchange boson, so Hc2{T) values may signifi- 
cantly exceed H^^^ as has been observed in other su- 
perconductors, for example PbMoeSs^i^. Thus, Eqs. 
(|8|)-(|9]) may capture the qualitative behavior of Hc2{T) as 
functions of materials parameters if they are expressed in 
terms of the observed quantities, and the Zeeman terms 
are renormalized by strong coupling effects. 



VI. DISCUSSION 

Our high-field data obtained on a Nd(Oi_a;Fa^)FeAs 
single crystals show extremely high Hc2 values, con- 
sistent with previous low-field measurements on sin- 
gle crystals^i^i^i^i^ and high-field measurements on 
polycrystals^. Our values of the temperature- 
dependent Hc2 slopes, H'JI ^ S - 10 T/K and H'^^ ^ 
1.2 — 5 T/K measured from the resistivity onset of su- 
perconductivity are considerably higher than the val- 
ues of Hc2 extracted from calorimetry measurements on 
Nd(Oi_a;Fa;)FeAs single crystals^^, which may be af- 
fected by strong paramagnetism of Nd^+ ions at high 
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FIG. 16: Fits of Eqs. {|8|-(|9| to the experimental data shown in Figs. [8] and [9] for different pairing scenarios: w > 0, 
All = A22 = 0.5, A12 = A21 = 0.25, r]± - D^^^~^ /Df^ = 0.08, ri\\ = [D^^^~^ D^^^'^ / d[^^ d[''~^]^^'^ = 1, and neghgible paramagnetic 
effects (a); ^ < 0, An = A22 = 0.49, A12 - A21 = 0.5, ry^ = i:>^"V^f ^ = 0-007, r]\\ = [D^^""^ D^^"'^ / d[^^ d["^Y^'^ = 1, and neghgible 
paramagnetic effects (b); 11; > 0, An = A22 = 0.5, A12 = A21 = 0.25, ry± = D^^V^f^ = 0.08, ri\\ = [D^^^~^ D^^^'^ / d[^^ d[^'^Y^'^ = 1, 
= d[^^/Dq = 10, and d\\ = [D^^^^D^^'^ /D^^^D^^^Y^^/Dq = 0.7 (c). 



field^. We should also mention several issues, which 
may complicate accurate extraction of the anisotropy 
in Nd(Oi_a;Fa;)FeAs single crystals grown by the flux 
method, which usually results in intergrowths and stack- 
ing faults parallel to the ab plane. Such planar defects 
may increase the apparent anisotropy parameter 7 and 
give rise to the local modulations of order parameter 
along the c-axis. The latter may broaden the resistive 
transition at Tc and the upward curvature of Hc2 {T) near 
Tc, which could mask the two-band effects. 

Because of very high values of H'JI and our high 
field measurements up to 60 T are still mostly limited to 
the temperatures not too far from Tc, particularly for the 
field parallel to the ab plane. As a result, it is hard to un- 
ambiguously evaluate the values i^c2(0), as Figs. [T6lf a-c) 
illustrate. Yet we can evaluate the in-plane and the c-axis 
GL coherence lengths, = (^o/27ri^^2 ^c)^/^ - 2.3 nm 
and = iah = 0-26 nm for 7 = 9, i^^2± = 1-2 T/K 
and Tc = 50 K. These values of and ^c are not too 
different from those for YBa2Cu307_a;. The fits shown 
in Figs. [T6lf a-c) also indicate the importance of paramag- 
netic effects in i^c2 (T) , particularly the fact that extrapo- 
lations of i^c2 to low temperatures based on the measured 
slope i^^2 and the WHH formula Hc2 ^ O.lTcH'^^^ may 
grossly overestimate Hc2{0). At the same time, the im- 
portance of paramagnetic effects in oxypnictides might 
open up possibilities of observing first order phase tran- 
sitions or perhaps the Fulde-Ferrel-Larkin-Ovchinnikov 
inhomogeneous states at very high magnetic fields. 

The temperature dependence of H^{T) shown in Fig. [9] 
is rather different from the conventional one-band WHH 
behavior, which may be characteristic of many oxypnic- 
tides. This confirms our earlier qualitative conclusion 
inferred from indirect measurements of H^{T) on poly- 
crystalline oxypnictides^. One possibility to explain the 



significant upward curvature of H^{T) is based on the 
two-band model. This conclusion is also consistent with 
the temperature dependent mass anisotropy parameter, 
7(T), reminiscent of the behavior of MgB^^i^. Yet 
many important points remains controversial. First, to 
explain the observed upward curvature of i^^(T), we had 
to assume a rather high difference in the intraband dif- 
fusivities, D2 ~ (0.1 — 0.01)T)i (depending on the pair- 
ing scenario), considerably higher than the effective mass 
difference for F and M electron pockets predicted by ab- 
initio calculation^iiiii^ii^. One has therefore to assume 
that scattering on impurities or strong magnetic excita- 
tions may produce such big differences in Di and D2. 

Our data enable us to make further qualitative con- 
clusions regarding the s-wave two-band pairing scenar- 
ios with and without the interband tt shift. The first 
one is less susceptible to the pairbreaking effect of in- 
terband impurity scattering, as follows from Eqs. ([5]) 
and (|6|). Indeed, the fact that our Nd(0,F)FeAs single 
crystal has the same Tc as sintered multiphase polycrys- 
tals, seems to suggest that the impurity scattering does 
not suppress Tc much. The latter may also indicate that 
oxypnictides are effectively in the clean limit just because 
of their short coherence lengths, < ^5 where i is the 
mean free path. Yet pairing is certainly more suscep- 
tible to the pairbreaking interband impurity scattering. 
Another indication that the s-wave pairing without the 
interband tt shift appears to fit our data better comes 
from the fact that, in order to explain the observed up- 
ward curvature of H^(T) in Fig. [T6lfb) in the model, 
we have to assume the ratio Di / D2 to be about an order 
of magnitude smaller that D2/D1 ~ 0.1 required for the 
more traditional model of strong intraband pairing and 
weaker interband pairing. These constraints may reduce 
the number of possible scenarios of superconductivity in 
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the oxypnictides. 

As far as vortex dynamics in Nb(F,0)FeAs is con- 
cerned, our data presented in Figs. [3] and [H show that 
high magnetic fields parahel to the c-axis mostly shift 
the resistive transition without significant broadening 
of the pxxiT) curves, while H parallel to the ab plane 
changes the shape of pxx{T)^ but certainly not as strong 
as for high- Tc cuprates^^. This suggests that despite the 
rather high values of 7, thermal fiuctuation of vortices 
in Nb(F,0)FeAs single crystals are weaker than in most 
anisotropic layered cuprates, like Bi-2212. This conclu- 
sion is consistent with the data of other groups^^ and 
with the relatively high activation energy £^0 3 x 10^ K 
extracted from our transport measurements. The effect 
of thermal fiuctuations is also quantified by the Ginzburg 
parameter, Gi = (27r/cBT'c/ioAo/0o^^c)^/2, where Aq is the 
in-plane London penetration depth. Taking Aq = 200 
nm, = 0.26 nm and Tc = 49 K, we obtain Gi 10~^, 
of the order of the typical Ginzburg number for YBCO. 



For two-band superconductors, the above estimates of 
Gi remain qualitatively the same if 7 and ^ are taken for 
the band with the minimum effective mass or maximum 
electron mobility^^. By contrast, Gi in La(0,F)FeAs is 
close to Gi in MgB2, about 30 times smaller than Gi in 
Nd(0,F)FeA*^^. 

In conclusion, our high-field magneto-transport mea- 
surements on single crystal Nd(Oo.7oFo.3o)FeAs have re- 
vealed very high upper critical fields and their anoma- 
lous temperature dependencies. We also observe a pro- 
nounced thermally-activated fiux fiow resistivity which 
may indicate rich vortex dynamics in single layer oxyp- 
nictides. 
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